Ferrofluids as thermal ratchets 
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Colloidal suspensions of ferromagnetic nano-particles, so-called ferrofluids, are shown to be suit- 
able systems to demonstrate and investigate thermal ratchet behavior: By rectifying thermal fluc- 
tuations, angular momentum is transferred to a resting ferrofluid from an oscillating magnetic field 
without net rotating component. Via viscous coupling the noise driven rotation of the microscopic 
ferromagnetic grains is transmitted to the carrier liquid to yield a macroscopic torque. For a simple 
setup we analyze the rotation of the ferrofluid theoretically and show that the results are compatible 
with the outcome of a simple demonstration experiment. 

PACS numbers: 5.40.-a, 82.70.-y, 75.50.Mm 



To extract directed motion from random fluctuations 
is a problem at the heart of statistical mechanics with 
long history 0, 0, 0, • On the one hand it is crucial 
for basing the second law of thermodynamics on statisti- 
cal reasoning and for depriving various Maxwell's demon 
type devices of their mystique. On the other hand it 
is closely linked to practical issues like the efficiency of 
heat engines. The problem has gained renewed atten- 
tion [E H> due to its possible relevance for biologi- 
cal transport 0, and the prospects of nano-technology 

A particularly clear example of a so-called ratchet is 
provided by an over-damped particle in a time depen- 
dent on-off spatially sawtooth potential subject to some 
random noise 13] . In the stationary (on or off) state, 
thermodynamic equilibrium prevails and detailed bal- 
ance prohibits directed transport. Complementary, with 
the potential switched on and off regularly diffusion and 
relaxation combine to yield a non-zero average particle 
drift. Other means may be used to drive the system 
away from equilibrium. In any case a noise driven drift 
is to be expected generically under non-equilibrium con- 
ditions where, however, the specification of the value and 
even the sign of this drift can be rather subtle . 

In the present letter we show that colloidal suspensions 
of ferromagnetic nano-particles, so-called ferrofluids [l^ . 
are ideal systems to test theoretical predictions on fluc- 
tuation driven transport experimentally. This is due to 
three main reasons: First, the small size of the ferro- 
magnetic grains (~ lOnm) implies that their dynamics is 
strongly influenced by thermal fluctuations [l5j • Second, 
spatially periodic time-dependent potentials for the ori- 
entation of the particles can be easily realized by external 
magnetic fields. Third, directed orientational transport 
manifests itself as systematic rotation of the ferromag- 
netic particles which can be easily detected from the re- 
sulting macroscopic torque on the carrier liquid. We show 
that our theoretical consideration are in agreement with 



a simple demonstration experiment. 

Ferrofluids combine the hydrodynamic properties of 
Newtonian fluids with the magnetic behavior of super- 
paramagnets |l4j . Many of their fascinating properties 
stem from the viscous coupling between the rotation of 
the magnetic grains and the vorticity of the hydrody- 
namic flow. A direct way to set this coupling into action 
is with a rotating external magnetic field as used, e.g., 
to spin up a ferromagnetic drop floating under hydraulic 
zero-gravity [Til Il7| . Complementary, a rotational fer- 
rofluid flow exposed to a static magnetic field exhibits an 
enhanced shear viscosity 0]. Various unusual effects 
such as "negative" rotational viscosity 19], magneto- 
vortical resonance [2fj|. and anomalously enhanced ac- 
response due to coherent particle rotation |2J rely on the 
exchange of angular momentum between rotating parti- 
cles and an oscillating magnetic field. In these cases the 
involved non-zero flow vorticity is crucial for explicitly 
breaking the symmetry between clockwise and counter- 
clockwise particle rotation. 

In contrast to these situations we demonstrate the 
transfer of angular momentum from an oscillating mag- 
netic field to a ferrofluid at rest. A crucial ingredient of 
the operating mechanism is a sufficient impact of thermal 
noise due to random collisions between the ferromagnetic 
grains and the molecules of the carrier liquid. 

We consider a ferrofluid in a horizontal, spatially ho- 
mogeneous magnetic field H composed of a constant part 
H x parallel to the x-axis and an oscillatory part H y f(t) 
along the y-direction. The precise form of the time mod- 
ulation f(t) will be fixed below. As common for diluted 
ferrofluids 14, we model the ferromagnetic particles 
as non-interacting spherical equal sized dipoles of volume 
V and magnetic moment m [2^. The external magnetic 
field gives rise to a potential U = — m • H for the di- 
rection of the magnetic moment. This direction in turn 
is assumed to be tightly co-aligned with the orientation 
of the magnetic particle. A reorientation of the mag- 
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netic moment hence requires a rotation of the particle 
against the viscosity r\ of the carrier liquid characterized 
by the Brownian relaxation time Tg = ZriV/ksT with T 
denoting the temperature of the liquid. We use the di- 
mensionless Langevin parameters a x ,y = mH XiV /kBT to 
measure the magnetic field strengths and scale time by 
twice the Brownian relaxation time. As time dependence 
in the ac part of the field we choose 



f{t) = cos(wi) + as\a(2ujt + f3). 



(1) 



The parameters a x , a y , a, u>, and (3 are easily controlled 
in experiments. 

Parameterizing the orientation of a magnetic dipole by 
two angles, e = m/m = (sin(9cos0,sin#sin</>, cos#) the 
potential U generated by the above specified magnetic 
field takes the form 

U(6,4>,t) = — sinO(a x cosfi + a y f(t) sin^). (2) 

Neglecting thermal fluctuations, the particle orientation 
is governed by an overdamped relaxation dynamics in 
the force field —VU. Accordingly, the orientation e(t) 
tends to align with the magnetic field. Since the field 
is restricted to the x-y-plane, 9(t) — > tt/2 for t — > oo. 
Within the plane, the angle between the field and the x- 
axis is always less than 7r/2, implying the same for cj)(t) 
when t — > oo, see Fig^ i.e. no average rotation or torque 
can arise. 

Things change in the presence of thermal fluctuations 
causing stochastic transitions between the deterministic 
solutions. As elucidated qualitatively with Fig^the dy- 
namical asymmetry induced by a x ^ and a ^ gives 
rise to slightly different probabilities for a 27r-phase-slip 
in the "forward-^ direction" and "backward-;/) direction" 
respectively. This in turn results in an average rotation of 
the particle, a manifestation of the ratchet effect [si Irl E^l] 
for the orientational motion of the ferromagnetic grains. 
Similar to situations with additive driving |24ll25l l26. 27] 
the spatial symmetry of our potential U (9, <j>, t) requires 
a sufficiently complex time modulation function f(t) for 
noise induced transport to occur. 

A more quantitative analysis of the effect can be made 
on the basis of the Fokker-Planck equation for the prob- 
ability distribution P(9, <j> : t) of the particle orientation. 
It is given by (see e.g. 28] ) 



d t P = V(P V U) + V 2 P, 



(3) 



where V denotes the angular part of the Nabla-operator 
in spherical coordinates. 

We have solved this equation numerically by expand- 
ing P(9,4>,t) in spherical harmonics. From the solu- 
tion we can determine the average orientation (e) = 
J d<j)d8 sin 6 eP(9, <j), t) and the average torque N p = 
m(e) x H acting upon an individual particle. For symme- 
try reasons only the z-component of the torque is differ- 
ent from zero. Due to the viscous coupling between the 




FIG. 1: Space-time plot of the potential U(8 = n/2,(j),t) 
from Eq.(|5J for a x — 0.3, a y = 1, a = 1, ui = 1, and f3 = 0. 
Blue and red regions correspond to small and large values of U 
respectively. In the long-time limit, the deterministic dynam- 
ics approaches 9(t) = tv/2 and a periodic (f)(t), represented by 
either of the full black lines. In the presence of thermal noise, 
transitions between these deterministic solutions become pos- 
sible schematically indicated by the dashed lines. The spatial 
asymmetry and temporal anharmonicity of the potential re- 
sults in slightly different rates for noise induced increments 
and decrements of cf> respectively. As a result a noise driven 
rotation of the particles arises. 



ferromagnetic grains and the carrier liquid the individual 
torques combine to yield a macroscopic torque per fluid 
volume 



N = nN p = fi Q M s (e) x H = fi Q M x H. 



(4) 



Here n denotes the number density of the particles, M = 
M s (e) the magnetization of the ferrofluid, M s = mn/fio 
its saturation value and the permeability of free space. 
A non-zero torque after averaging (@J over one period of 
the external ac-driving is the macroscopic signature of a 
directed orientational transport. 

An explicit expression for the torque can be obtained 
from an approximate solution of the Fokker-Planck equa- 
tion © by adopting the effective field method |2^. It 
gives rise to a relaxation equation |3(| for the average 
orientation (e) which for a fluid at rest reduces to 



(1 )0ij 

OLcS 



( 3(f)_ 1} ^M 
a e ff (e) (e) 



(a e ff— a)j 
(5) 



where (e) denotes |(e)| and the so-called effective field 
c* : ff((e)) is given by the inverse of the Langevin function 
(e) = C(a) = cotha — 1/a. The numerical integration 
of Eq.© yields rather accurate approximations for the 
orientation (ei)(t), which deviate from the numerical so- 
lution of the Fokker-Planck equation by less than a 
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FIG. 2: Sketch of the experimental setup. ( A side view, 
B top view). A hollow plastic sphere (1) is filled with fer- 
rofluid and suspended on a thin Kevlar fiber. The time depen- 
dent magnetic field in {/-direction is generated with a pair of 
Helmholtz coils (2), the static field in ^-direction stems from 
a commercial electro-magnet (3). Noise assisted transfer of 
angular momentum from the magnetic field to the ferrofluid 
manifests itself in a rotation of the sphere. 



few percent. Nevertheless the final values for the time 
averaged magnetic torque 0} differ by a factor between 
2 and 3. This is due to the fact that the time average 
N z is much smaller than the typical values of the time 
dependent component N z of the torque. 

For sufficiently weak magnetic field we may use C{a) ~ 
a/3 which makes eq. (J5J linear in (e). It is then eas- 
ily solved for a time dependence of the field a as spec- 
ified by eq.£[J. The resulting time averaged torque |Q}, 
however, vanishes identically. Exploring the non-linear 
regime pcrturbatively by using L(a) ~ a/3 — cv 3 /45 we 
find to leading order for the time-averaged ^-component 
of the torque 
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a u} 2 [u> cos (3 + 2 sin j3) 
30 (l + w 2 )(4 + u; 2 ) 2 ' 



(6) 



where x denotes the magnetic susceptibility of the fer- 
rofluid. This expression provides a useful approximate 
formula for the nontrivial dependencies of the torque on 
the parameters of the problem. In particular it shows 
that both the static magnetic field component a x and 
the anharmonic time dependence of the oscillatory 
component are essential for directed rotational transport 
to occur. The effective field method may be an useful 
tool for the analysis of other ratchet systems as well [3l| • 
To verify our theoretical findings we have designed a 
simple demonstration experiment as sketched in Fig|21 
The setup is a torsion balance similar to those used in 



string galvanometers. A hollow plastic sphere with in- 
ner diameter 16 mm was filled with a ferrofluid APG 
933 (Ferro-Tec), with the following specifications: den- 
sity p = 1, 100 kg/m 3 , susceptibility \ — 1-09, satu- 
ration magnetization M s = 18 kA/m, dynamic viscos- 
ity 77 = 0.1 Pas. The container was suspended on a 
Kevlar fiber with a length of 20 cm and 10 fim diameter. 
The oscillatory magnetic field was generated by a pair 
of Helmholtz coils of 100 windings each via a computer 
generated signal of the form specified in Eq.Q. The 
RMS-field strength in the center of the coil amounted 
to about H y ~ 2.1 kA/m. The static field component 
H x ~ 9 kA/m was generated by a commercial electro- 
magnet (Bruker). The applied frequency was v = 200 Hz. 
In the choice of the frequency and the ratio between the 
amplitudes of the static and oscillating field we where 
guided by our numerical solutions of the Fokkcr-Planck 
equation ©. The dimcnsionless parameters describing 
the setup are a x — 1.72, a y — 0.4, a = 1. The dimen- 
sionless frequency ui depends on the Brownian relaxation 
time tb which was used as a fit parameter. 

The main qualitative result of the experiment is the 
unambiguous demonstration of the proposed thermal 
ratchet effect: After switching on the fields the ferrofluid 
sphere immediately starts to rotate. Switching off cither 
the static field a x or the modulation amplitude a the 
torque disappears. A reversal of the sign of either a x or 
a causes a reversal of the sample's rotation sense. 

A more quantitative comparison between theory and 
experiment was difficult for the following reason: Our 
setup was designed to be as sensitive as possible with re- 
spect to the appearance or not of a rotation of the drop. 
The actual rotation speeds of the sphere then became 
usually so large that the field had to be switched off af- 
ter a short time in order not to destroy the experiment. 
We hence resorted to the following semi-quantitative pro- 
cedure to determine the torque: Starting with the sphere 
at rest and then letting act the magnetic field for a fixed 
amount of time (a few seconds) , the resulting final rota- 
tion speed of the sphere was used as a measure for the 
torque. In this way we have experimentally recorded the 
torque (in arbitrary units) as a function of the phase an- 
gle (3 in £[}. The results are shown in figCU together with 
a fit to our approximate theoretical result given in eq.JH}- 

Theory and experiment are compatible if the Brow- 
nian relaxation time is chosen as tb — 1.8 ms. Using 
t~b = 'ir/V/ksT this corresponds to a particle diameter 
of about 35 nm, which exceeds the typical size by about 
a factor of 3. However, in real ferrofluids size, shape, 
and magnetic moment of the ferromagnetic grains vary 
quite significantly (see e.g. Q, page 41), a situation 
only poorly described by our model involving a single re- 
laxation time. Moreover, since the transferred angular 
momentum increases with the particle diameter |3jJ the 
above value for tb is likely to characterize the largest 
grains in the population rather than the average. 
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FIG. 3: Magnetic torque transferred to the ferrofluid as func- 
tion of the phase angle f3 for a x = 1.72, ct y — 0.4, a = 1. 
The symbols are experimental results obtained with the setup 
shown in FigEI The solid line is a fit to the approximate an- 
alytical expression (|SJ, N z = A(uj cos ft + 2sin/3) with the 
amplitude A and the frequency u> as fit parameters. The 
value uj ~ 4.41 obtained translates into a fit for the Brownian 
relaxation time of tb — 1.8 ms. 

In conclusion we have shown that by rectifying rota- 
tional Brownian motion angular momentum can be trans- 
ferred from an oscillating magnetic field to a ferrofluid 
at rest. A unique feature ferrofluids offer in comparison 
with other experimental realizations of ratchets is the 
combined action of many individual nano-scalc ratchets 
to yield a macroscopic thermal noise induced transport 
effect. More precise experiments are needed to verify 
quantitatively our theoretical predictions about the de- 
pendence of the torque on the magnetic field strength 
and the frequency of the external driving. 

Many other investigations using ferrofluids as thermal 
ratchets are conceivable. A particularly interesting line 
is linked with the occurrence of inversion points Pi of the 
rotation, see Fig. |31 When keeping j3 fixed to such an 
inversion point and varying instead another parameter 
of the system the rotational speed of the magnetic grains 
will generically exhibit a change of sign as a function of 
this parameter as well. Since a x ,y or the dimensionlcss 
lo (due to its dependence on tb) depend on the size of 
the magnetic grains in a polydisperse ferrofluid under 
the same experimental conditions the larger and smaller 
particles may rotate in opposite directions. 
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Blunder, I. Rehberg, and M. Raible are gratefully ac- 
knowleged. This work was supported by Deutsche 
Forschungsgemeinschaft under SFB 613, SFB 277, FOR 
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